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Abstract 

We present the results of first-principles cluster calculations of the electronic 
structure of La2Cu04. Several clusters containing up to nine copper atoms 
embedded in a background potential were investigated. Spin-polarized cal- 
culations were performed both at the Hartree-Fock level and with density 
functional methods with generalized gradient corrections to the local density 
approximation. The distinct results for the electronic structure obtained with 
these two methods are discussed. The dependence of the electric-field gradi- 
ents at the Cu and the O sites on the cluster size is studied and the results are 
compared to experiments. The magnetic hyperfine coupling parameters are 
carefully examined. Special attention is given to a quantitative determination 
of on-site and transferred hyperfine fields. We provide a detailed analysis that 
compares the hyperfine fields obtained for various cluster sizes with results 
from additional calculations of spin states with different multiplicities. From 
this we conclude that hyperfine couplings are mainly transferred from nearest 
neighbor Cu 2+ ions and that contributions from further distant neighbors are 
marginal. The mechanisms giving rise to transfer of spin density are worked 
out. Assuming conventional values for the spin-orbit coupling, the total cal- 
culated hyperfine interaction parameters are compared to informations from 
experiments. 
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I. INTRODUCTION 



A large amount of information about both static and dynamic properties of materials that 
exhibit high-temperature superconductivity, is available through experimental techniques 
like nuclear quadrupole and magnetic resonances (NQR and NMR). In particular, the electric 
field gradients have been determined for a variety of nuclei with great accuracy! These 
quantities are ground-state properties of the solid and depend sensitively on how charge 
is shared among the atomic and molecular orbitals. Measurements of the Knight shift 
and the various nuclear elaxation rate tensors as a function of temperature, doping, and 
orientation of the applied field with respect to the crystalline axes, provide insight into the 
static spin density distribution and the low frequency spin dynamics of the electrons in the 
normal as well as in the superconducting stateB. Owing to the abundances of nuclei with 
suitable nuclear magnetic and quadrupole moments, most of these quantites can be studied 
on different sites in the unit cell. This allows one to distinguish static and dynamic features 
in the Cu02 planes from those in the interlayer region. 

The analysis of numerous early NQR and NMR experiments focussed on the question 
whether a one-component model is sufficient to describe the low-energy excitations of the 
quantum fluid in the high-temperature superconductors or whether a two-component de- 
scription is needed. Knight shift and relaxation rate data for the Cu sites in YBa2Cu307 
that were analyzed on the basis of the hyperfine coupling of a Cu 2+ ion in an axial sym- 
metric surroundig, seemed to be controversial. Mila and RicJ proposed then that besides 
the usual on-site hyperfine interactions at a Cu nucleus an additional hyperfine field should 
be considered which is transferred from neighboring Cu ions. They gave a consistent ex- 
planation of the unusual combination of anisotropies of the Cu Knight shifts and relaxation 
rates in YB^CusOy within a one-component theory. An estimation of the strength of the 
various hyperfine parameters with the use of a quantum chemical model led to good agree- 
ment with the values deduced from experiment. Monien, Pines, and Slichteii analyzed a 
more complete set of data from which they inferred the relevant hyperfine couplings and the 
associated fluctuation times. The ensuing constraints on theoretical one-component models 
for the planar magnetic excitations comprise a large transferred hyperfine field between ad- 
jacent Cu 2+ spins, which, moreover, must be very nearly aligned antiferromagnetically, and, 
in addition, an almost perfect cancellation between the parallel component of the dipolar 
hyperfine tensor and the total transferred field (A\\ + 4B w 0). Provided that the spin sus- 
ceptibility is isotropic, it is only with the latter condition that the vanishing of the parallel 
component of the spin contributions to the Knight shift can be explained. 

A remaining intricate problem whose explanation seemed to require a two-band model, 
was the relaxation behavior of the planar oxygen nuclei which was strikingly different from 
that of the coppers. Shastryi, however, pointed out that spin density is also transferred to 
the planar oxygens. A transferred hyperfine field couples a nucleus to spin fluctuations on 
neighboring Cu sites. This implies that the spin-lattice relaxation time T\ depends on form 
factors given by the local geometry of lattice sites of the nucleus and its neighbors. Since 
the form factor entering l7 Ti vanishes at the commensurate antiferromagnetic wave vector 
(it /a, 71 /a) where that for 63 Ti reaches a maximum, the variant temperature dependence of 
17 Ti and 63 Ti could be understood from the presence of strong antiferromagnetic correlations 
between the copper spins in the planes. Thus a coherent picture of the low-energy spin 
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dynamics in the Cu02 planes seemed to be established^ and numerous NMR and NQR 
experiments were analyzed on the basis of these hyperfine spin Hamiltonians and a one-band 
model. Pines and coworkers developedB a phenomenological expression for the dynamic spin 
susceptibility and provided!'! a quantitative description of a variety of NMR data. 

New doubts about the adequacy of a one-component description of the spin fluid arose 
when the results of neutron scattering data became available. NMR and neutrons both probe 
the spin susceptibility. In Lai.seSro.uCuCU, a considerable degree of incommensuration of 
the peaks in the spin-fluctuation spectrum was observed^. As a consequence, the NMR 
relaxation time l7 Ti is expected to exhibit an anomalous temperature dependence which is 
not seen experimentally. This and other contradictions between the results of NMR and 
neutron-scattering experiments led Zha, Barzykin, and Pines0 to a critical re-examination 
of the hyperfme spin Hamiltonian. They showed that it is possible to reconcile NMR and 
neutron-scattering experiments on both La 2 _ 3; Sr, E Cu04 and YB^CusOg+z within a one- 
component theory by introducing a transferred hyperfine coupling between the O nuclei 
and their next-nearest-neighbor Cu 2+ spins. In addition, the analysis exhibited that the 
hyperfine coupling transferred to the Cu nuclei changes as one goes from the L^-xSr^CuC^ 
to the YBa2Cus06+x system, and is moreover comparatively sensitive to hole doping in the 
former system. 

Despite the rich information of NMR and NQR results about the nature of the spin fluid 
and the low-energy quasiparticles, there exist only few theoretical first-principles approaches 
which address the determination of electric-field gradients and magnetic hyperfine interac- 
tions. For the YBa2Cu306+x system, electric-field gradients (EFGs) at the various nuclear 
sites have been obtained by Das and coworkers^TEJ with ab-initio cluster calculations us- 
ing the unrestricted Hartree-Fock (UHF) method, and by Schwarz and coworkersElJll who 
employed the full-potential linear augmented-plane-wave method within the local density 
approximation (LDA). Results of Hartree-Fock (HF) calculations have also been published 
by Winter^. Apart from one exception, the EFGs calculated with these three different 
methods more or less agree and reproduce the experimental data quite satisfactorily. The 
exceptional case is the EFG at the planar Cu sites. Using a large cluster that contained 74 
atoms, we recently performed! calculations with the density functional (DF) method and 
obtained values for the Cu EFGs which are in better agreement with the experiments. These 
calculations were non-spin polarized and no information on magnetic hyperfine interactions 
is available. The energy difference between the singlet and triplet states of two coppers in 
adjacent Cu0 2 planes was studied^ in a small cluster containing 2 Cu, 4 Y, and 8 O atoms. 

For the La2Cu04 system, UHF calculations have been reported by Sulaiman et al.0 who 
obtained EFG values at various nuclei as well as the contact and dipolar hyperfine fields at 
Cu. Martin and Hay0 presented HF results of the electronic structure of Cu06 clusters in 
neutral, p-doped, and n-doped states and studied the influence of correlation effects using 
the technique of configuration interactions. Martini! calculated HF values for EFGs in 
Cu0 6 and CuoOn clusters and investigated the change in the NQR spectra upon doping. 
A comparison^ of HF and LDA calculations for a CuOg cluster showed large differences in 
the EFG values and the hyperfine coupling parameters. 

In this paper, the results of extensive cluster studies for the pure La2Cu04 system are 
presented. Spin-polarized calculations at the HF level and with the DF method with gradient 
corrections to the correlation functionals have been performed for clusters comprising n = 
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1, 2, 4, 5, and 9 Cu atoms in a plane. The resulting electronic structures, the charge and 
spin distributions, EFGs and magnetic hyperfine interactions have been analyzed. In Sec. ||, 
details on the clusters, their embedding in a lattice of point charges, the basis sets and the 
theoretical methods are given. Results for the CuOg 0- ion are presented in Sec. [Hj and 
compared to the predictions for the hyperfine couplings of a Cu 2+ ion in a single-electron 
picture. The distinct results of HF and DF methods are discussed in detail. Section [IV] is 
devoted to the investigation of the EFG values for Cu and the planar and apical oxygens and 
their dependence on the cluster size. The results for the magnetic hyperfine interactions are 
given in Sec. [V|. Particular emphasis is put on investigations of contributions to the hyperfine 
fields transferred from nearest and further distant copper ions. The origins of the various 
hyperfine couplings are worked out and explained in terms of localized atomic orbitals and 
molecular orbitals. The general mechanisms giving rise to spin density transfer are discussed. 
Sec. [VI| contains a summary and conclusions. 



II. CLUSTERS AND BASIS SETS 

The general idea of the cluster approach to electronic structure calculations is that the 
parameters that characterize a small cluster should be transferable to the solid and largely 
determine its properties. The essential contributions to EFGs and to magnetic hyperfine 
fields are given by rather localized interactions and therefore it is expected that these local 
properties can be determined and understood with clusters calculations. Approximations 
must be made concerning the treatment of the background that is employed to embed the 
cluster. Using as large a cluster as is possible is of course advantageous. It is necessary, 
however, that the results obtained should be checked with respect to their dependence on 
the cluster size. 

The clusters used in this work comprise copper atoms of the Cu02 plane together with 
the surrounding planar and apical oxygens. They were embedded in a lattice of point 
charges with charges +2, —2, and +3 according to the ionic charges of Cu 2+ , O 2- , and La 3+ , 
respectively. However, it has been pointed out in Ref.0 that bare positive point charges close 
to the border of the cluster present a too strong attraction for the electrons. It is therefore 
essential to replace them by appropriate pseudopotentials which guarantee the orthogonality 
of the diffuse electron wave-functions with the ion cores. 

The smallest cluster investigated is shown in Fig. |I|. All-electron basis sets were put on 
one Cu and six O atoms and pseudopotentials were used to represent the four Cu 2+ and 
ten La 3+ sites closest to the cluster border. In addition, these 21 atoms were surrounded 
by more than 2000 point charges. Some charges at sites far from the cluster center were 
adjusted such that the total charge of the system was zero. Next, the positions of some 
point charges were changed in such a way that the correct Madelung potential in the central 
region of the cluster was reproduced. 

All atomic positions were located according to the tetragonal phase with lattice con- 
stantsii a = 3.770 A and c = 13.18 A and with a Cu-O(a) distance of 2.40 A. 

For the cluster atoms the standard 6-311 G basis sets were employed. For Cu, this 
corresponds to the basis set developed by Wachters@, while for O it was given in RefH 
In Table [I] an overview is presented of the clusters used in this work with a listing of the 
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numbers of Cu atoms, of cluster atoms, of electrons, of basis functions, and of primitive 
Gaussian functions. 

The majority of the HF and DF calculations were spin-polarized. They were per- 
formed with the Gaussian94 and Gaussian98 programs0. The Vosko-Wilk-Nussair func- 
tional (VWN)0 was used as local density approximation (LDA). For the gradient corrections 



ized gradient approx- 
on YBa2Cu30y, we 



to the exchange and correlation functionals several forms for the genera 
imations (GGA) are available in the literature. As in our invest igat ions 
mainly used the formula proposed by Beckei§ in combination with the functional of Lee, 
Yang, and Pari@ (BLYP). For the smallest cluster, other functionals have been used as well, 
such as the X Q -LYP (XALYP) and the form proposed by Perdew and WangEH (BPW91). 
The individual integrals over the atomic orbitals for the electric field gradient were taken 
from calculations using the GAMESS-US program!!. For the final analysis a special pro- 
gram was designed^ which employed these integrals and the electronic wave functions from 
the Gaussian calculations . 

Throughout this work we use atomic units for the electric-field gradients and magnetic 
hyperfme interaction parameters. To distinguish these from other notations, we will use 
lower case letters and write, e.g., the hyperfme spin Hamiltonian for a nuclear spin / and a 
single electron spin S as 

H = ^2 I a A af3 Sp = h'jjh'je la a a p S p . (1) 

a/3 a/3 

The tensor a a p has the dimension of a density. In the NMR literature, different and some- 
times misleading units are in use. We note that 2//g/a 3 3 = 125.2 kG, which connects to the 
units used in RefH 



III. EMBEDDED CU0 6 -ION 

To investigate the (CuOg) -10 ion embedded in an environment of the La2Cu04 system, 
the cluster Cu0 6 /Cu 4 Laio shown in Fig. [I] was used. The results for the EFG component 
V zz , the core polarization, the dipole tensor, and the atomic spin density p at the Cu site 
are collected in Table [H]. 

Our results for V zz and adip obtained at the HF level agree with previous predictions@H 
There exist, however, distinct differences between the HF results and those calculated within 
the DF/GGA scheme. In HF, the EFGs are about 40 % larger than those obtained from 
DF calculations. A similar discrepancy between EFG values calculated with HF and DF 
methods has also been observed^ in calculations of the EFG at the planar Cu site in the 
YBa 2 Cu30y system. 

The results of different calculations within the DF approach, performed by using other 
exchange and correlation functionals (SVWN, XALYP, BPW91), however, agree among each 
other reasonably well. Increasing the quality of the basis set by including polarization (p) 
and diffuse (d) functions also does not alter the results substantially. 

The experimental valued of the quadrupole frequency at the 63 Cu is 33.0 MHz. This 
corresponds to a zz = 1.331 assuming a quadrupole moment Q( 63 Cu) = —.211 b. A more 
detailed comparison between experimental and theoretical values for the EFG is given in 



5 



Sec. |V[ The results for the core polarization and the dipole tensor will be discussed and 
compared to experiments in Sec. |V[ 

It is instructive to compare these ab-initio results for the embedded (CuOg) -10 ion with 
the predictions for a single-electron (or rather single-hole) picture of the Cu 2+ ion which 
have been given by Bleaney et al.S. If the d-shell is full except for a singly occupied ?>d x 2_ y 2 
orbital, the dipolar hyperfine field is given by 

alp = ~< r~ 3 > (2) 
and the core polarization is determined by a parameter k: 

a cp = -k < r~ 3 > . (3) 

The spin-orbit coupling will be discussed later. 

Inserting the DF-value ajL, = —3.526 from Table [IT] into Eq. (1) gives < r~ 3 > = 6.171. 
Using this number and a cp = —1.784 in Eq. (2) leads to k = 0.289. The corresponding 
HF-values yield < r~ 3 > = 7.737 and k = 0.455. 

Compared with the density functional ab-initio results for the embedded (Cu0 6 )~ 10 ion, 
the single-hole picture already gives reasonable results for adip and a cp . It fails, however, in its 
prediction of the EFG with a value of V zz = —4/7 < r~ 3 >. This discrepancy has sometimes 
been used to assign a fractional hole number to the 3d x 2_ y 2 orbital. The underlying physical 
picture, however, is completely misleading since the contributions to the EFG originate from 



various Cu shells as can be seen in Table |TJ. The cluster EFG value of 1.4 as obtained with 
the DF method, results from cancellations between relatively large individual contributions. 
This shows that the theoretical determination of EFGs is quite delicate. It is necessary to 
describe all electron shells of the atom accurately. 

At this point we would like to discuss the differences in the results as obtained in HF- ver- 
sus DF- theory by analyzing the details of the electronic structure. For each spin projection 
there exists a total of 23 molecular orbitals (MO) that can be formed as linear combinations 
of the Cu 3d and the six O 2p atomic orbitals (AO). Both HF and DF predict that the 
MO with highest energy is the anti-bonding hybridization between Cu 3d x 2_ y 2 and O 2p x 
and 2p y AOs. It is occupied for one spin projection (spin up) only and, in the following, 
energies will be measured with respect to the energy of this highest occupied molecular or- 
bital (HOMO). (The HOMO calculated with DF for a cluster comprising 9 Cu atoms will be 
shown in Fig. [5].) The corresponding spin-down orbital is the lowest unoccupied molecular 
orbital (LUMO) and lies at an energy which is at 1.4 eV for DF and even higher for HF. 
The energies of all 23 MOs are shown in Fig. |2| together with the contributions from the 
individual AOs. The length of the bars marks the squared values of the expansion coefficient 
of the MO with respect to the AOs whereby spin-up (spin-down) orbitals are denoted by 
solid (dotted) bars. 

The orbitals with second highest energy are anti-bonding linear combinations of the AOs 
3d 3z 2_ r 2 of the Cu and O 2p z of the apical oxygens. In the DF calculation, these are followed 
at energies around -1.5 eV by MOs that can be formed as anti-bonding hybridizations with 
the three other Cu 3d AOs. Deeper in energy, between -2.6 and -4.4 eV, are the MOs that 
are composed of oxygens only without contributions from Cu. They comprise the 2p z AOs 
of the planar oxygens and the 2p x and 2p y of the apex oxygens. At the HF level, these are 
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at energies (-1.1 to -1.8 eV) that are above those of the MOs with contributions from Cu 

xy,yz,zx- 

The MO with lowest energy in Fig. ^| is at -5.8 eV for DF. It is formed by the bonding 
hybridization between 3d x 2_ y 2 and O 2p x and 2p y . Further down in energy lie 2 x 21 orbitals 
formed by inner-shell AOs. Their energies and wave functions are only slightly different for 
the two spin projections though it is precisely this difference that will determine the core 
polarization and also contributes to the electric field gradient. Therefore it is essential to 
describe also inner-shell electrons with basis sets of sufficiently high quality 

The HF result suggests an almost purely ionic bonding with small overlap between Cu 
and O. In contrast, the DF results emphasize the covalent character of the bonds and an 
appreciable overlap. This varient description of bonding is also reflected in the width of the 
"d-band" (10 eV vs. 6 eV) and the localization of the atomic spin density at the Cu (0.90 
vs. 0.67). 

The contributions of the DF MOs shown in Fig. to the density of states is represented 
in Fig. |3|. The individual states were broadened by folding with a Gaussian function of 
half-width 0.21 eV. A comparison to the density of states obtained from band structure 
calculations0S shows an overall agreement even for the small cluster under consideration. 

The total occupation of the DF Cu atomic orbitals amounts to 8.92 which results from 
the following contributions: .08, .89, .96, .96, .92, .97 from 4s, 3d 3z 2_ r 2, 3d zx , 3d xy , 3d x 2_ y 2, 
3d xy with spin up and .07, .66, .96, .96, .27, 1.00 with spin down. The small population of 
the Cu 4s AO is due to the fact that the hybridization between the 3d3 Z 2_ r 2 and the 2p z of 
the apex oxygens also involves a minor admixture of 4s orbitals. The populations obtained 
from the analysis of the HF results are somewhat higher, in agreement with Sulaiman et 
alS 

In conclusion, we think that the Cu-0 bonding is better described by DF methods, which 
include part of the correlation, than with HF. This will gain further evidence from the results 
obtained with larger clusters that will be reported in the following sections. 

IV. ELECTRIC FIELD GRADIENTS 
A. Larger clusters 

To investigate the convergence of the calculated EFGs and magnetic hyperfine properties 
with respect to cluster size, calculations for clusters comprising n Cu atoms (n = 1, 2, 4, 5, 9) 
have been performed using the HF method as well as DF with GGA (BLYP-functional). 
The multiplicity m of the spin states was chosen to be maximal (m = 2n + 1). Lower 
multiplicities were also investigated for selected clusters (m — for n = 2 and m = 2, 4 for 
n — 5). On all atoms belonging to the cluster, a 6-311G basis set was employed. In Fig. |] 
the largest cluster used is represented. It comprises 9 Cu and 42 O atoms in the plane and 
is surrounded by 12 Cu and 50 La ions that were treated with bare pseudopotentials. The 
corresponding HOMO as obtained with DF methods is shown in Fig. |5|. 

In this section, the resulting EFG values are discussed whereas the magnetic hyperfine 
properties will be analyzed in Sec. [V|. 
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B. Electric Field Gradients at the Cu Site 



The results of Hartree-Fock and density-functional calculations for the V zz component 
of the EFG at the Cu site are given in Table |TV|. Since the Cu sites in the two clusters 
with even number of Cu atoms are not situated in the center, a small asymmetry rj results. 
Again, the large discrepancy between HF and DF results is obvious. Our HF values for 
small clusters are in agreement with those obtained in Refs. 

The variations of the values for the EEGs with respect to cluster size and multiplicity 
are within reasonable bounds considering the subtle cancellations of contributions from the 
various shells as can be seen in Table [TTT| . 

The experimental value!! for the 63 Cu quadrupole frequency is 33.0 MHz. This corre- 
sponds to a value for V zz of 1.331 a.u. using a quadrupole moment of 63 Q = —0.211 b, but to 
1.560 a.u. with 63 Q = —0.18b. The former value for 63 Q is from an analysis of Sternheimerll 
and is the one cited in the current NQR tables while the latter was determined by Stein et 
al.0 from a HF cluster calculation for cuprite. In relation to the large amount of NQR data 
available on cuprate superconductors, the availability of reliable values for the quadrupole 
moments of Cu would be very desirable. 

In a comparison with experiments it should further be noted that our calculations were 
performed for clusters with atomic positions corresponding to the tetragonal phase. More- 
over, spin states with lower multiplicities better account for antiferromagnetic fluctuations 
than high-spin states. In this respect, our calculated values are by 17 % too low, even with 
m Q = —0.2116. Nevertheless, it can be concluded that a valuable estimate for V zz can al- 
ready be gained with DF calculations with the cluster CuOe/C^Laio. This is of relevance 
for studies of the changes in the quadruole frequency upon doping. 



C. Electric Field Gradients at the O Site 

Although the clusters used in this report were constructed to investigate mainly the 
magnetic hyperfine interactions at the Cu site, the results obtained for the EFGs at the 
oxygen sites are also of interest. For the planar oxygen O(p) they are listed in Table [V]. 
Again, we find reasonable agreement between the experimental values and our calculations 
with the DF method. 

For O(a), the apical oxygen, the EFG is axially symmetric and the component V zz turned 
out to have the values .26, .21, and .19 for DF (.28, .26, and .19 for HF) for clusters with 4, 
5, and 9 Cu atoms, respectively. The experimental valued is .22. 



V. CORE-POLARIZATION AND TRANSFERRED HYPERFINE FIELDS 

A. Isotropic Contributions 

The magnetic hyperfine interaction at a nucleus can be decomposed into an isotropic 
part D and a traceless dipolar part. The former is given by the difference between the spin 
densities at the nuclear site R: 
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D(R) = ^fel^(R)l 2 -E 

6 \ ni m ' 



^(R)l 2 ) (4) 



where the sums extend over all occupied MOs. This quantity is called Fermi contact con- 
tribution if the singly occupied MOs contain AOs of s-symmetry centered at R. It is called 
core-polarization if the MOs are doubly occupied but with differing contributions from the 
s-like AOs. In the present case, the highest occupied molecular orbital is always composed 
of a hybridization of predominantly atomic Cu 3d x 2_ y 2 and O 2p XjV character. Therefore we 
prefer the use of the notation core-polarization in the following. 

The values for D(Cu) vary strongly with cluster size and position of a particular Cu atom 
in larger clusters. In the Cug042/Cui 2 La4o cluster, e.g., D = 0.922 for the central Cu which 
is surrounded by 4 nearest neighbors (NN), but D = —0.364 at the four equivalent corner 
sites which have 2 NN and D = 0.321 at the four sites on the edges with 3 NN. The complete 
set of data is compiled in Table [VI]. In Fig. ||, the values for D(Cuj) obtained with clusters 



of different sizes are plotted vs. N, the number of nearest neighbor Cu atoms of Cuj. In all 
calculations, the same basis set 6-31 1G was used. It is obvious that in a first approximation 
the increase of D(Cu) is linear in the number of nearest neighbor Cu atoms for both DF and 
HF methods. It is suggestive to attribute this change to a transferred hyperfine term b and 
to decompose D(Cu) into 

D(Cu) = a + bN (5) 

with preliminary values: a ~ —1.75 and b « 0.69 for DF/BLYP and a « —3.52 and b ~ 0.17 
for HF. 

For those oxygen sites that are at the border of the cluster and that are therefore bounded 
to a single Cu atom, a value D(O) ~ 0.64 is found. For oxygens bridging two Cu atoms, 
however, we get D(O) m 1.24 with the HF values 0.60 and 1.11 being only slightly lower. 
Therefore, also the core polarization at O is determined by a transferred hyperfine field from 
its neighboring Cu atoms: 

D(O) = c N (6) 

with a preliminary value c ~ 0.6. 

At the apex oxygens, the spin density is marginal (on the order of 1 %) and the core 
polarization amounts to 0.05. 

In retrospective, some of the early problems in analyzing NQR and NMR data on cuprate 
superconductors can be traced back to the fact that the single-hole model developed in 
Ref.0 for an isolated Cu 2+ ion in a crystal field, cannot be applied to correlated Cu ions in 
a Cu02 plane. To be precise, it works reasonably well for local quantities like the dipolar 
coupling tensor and the spin-orbit interaction which will be discussed later. Their values are 
mainly determined by the HOMO and also do not depend too sensitively on the method of 
treating the many-electron system. Even HF calculations give results of the correct order of 
magnitude. The core-polarization, however, is not a local affair. As Fig. ^| exhibits, hyperfine 
fields transferred from the neighboring copper ions have a drastic effect on D(Cu). Its value is 
negative for an isolated copper ion but is increased by contributions from neighboring ions. 
At the HF level, there is a small hybrization and the transfer is small. DF calculations, 
however, render a stronger covalency and therefore the spin density is more delocalized 
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which implies a larger transfer. With contributions from four NN Cu ions, D(Cu) becomes 
positive as was first recognized by Mila and Rice!. 

The question then arises whether also contributions from further distant Cu atoms will 
contribute to the hyperfine properties at Cu or O. This is of particular relevance since to 
reconcile NMR and neutron-scattering experiments within a one-component spin model, 
Zha et al.0 have advocated a hyperfine field at oxygen transferred from second nearest Cu 
neighbors. This would lead to a reduction of the form factor 17 F(q) around the wave-vector 
q — (it /a, it /a) and to a modification of the spin-lattice relaxation rate 17 T 1 _1 . 

To investigate this problem of hyperfine fields transferred from further distant neighbors 
in detail, we performed also DF cluster calculations for spin states with lower multiplicity. 
In the cluster Cu5026/CusLa34 states with multiplicities m = 4 and m = 2 were obtained. 
For the former the spin density at the central Cu atom is of opposite sign to that of the four 
neighboring sites. For m = 2, it is negative at the central and one edge Cu atom. In Fig. [7] 
the corresponding differences in the spin densities along the Cu-O-Cu-O-Cu axes are shown. 
They peak at positions where the squares of the 3d a .2_ 3/ 2 orbitals have their maxima. The 
core-polarization then yields a change close to the Cu nucleus. In the "antiferromagnetic" 
cluster that is obtained with m = 4, the atomic spin density at the central Cu is negative 
but its magnitude is smaller than those at the four adjacent coppers (the detailed results 
are given in Table |V11| ). This is due to the finite size of the cluster. There is an excess of 
positive spin density in the central region which is also reflected by the atomic spin densities 
at the four oxygens that are bridging two antiferromagnetically coupled Cu atoms. There, 
a value of 0.026 is obtained instead of the cancellation expected for a large system. 

To account for these finite size effects, we have carefully analyzed all results given in 
Tables |V| and [VI I| . It turned out that the values of the core-polarization at the oxygens 
are proportional to the sum of atomic spin densities at the adjacent Cu ions 

D(O) = 7 E P(Cuj) (7) 

jeNN 

where the sum extends over one or two NN. This is illustrated in Fig. [8|. The negative values 
result from the spin states with lower multiplicities. The deviations from the linear behavior 
for some small D(0)-values is due to errors in forming the difference between large p(Cuj) 
values with opposite signs. 

Similarly, the core-polarization at the Cu ions can well be reproduced by decomposing 
it according to 

D(Cui) = ap(CuO +(3 E p(Cuj) . (8) 

jeNN 

In Fig. |] the various values Du(Cui)/p(Cu(i)) are plotted as a function of the sum of NN 
atomic spin densities divided by p(Cuj). 

With all these data available, we are now in a position to study the intrinsic and trans- 
ferred hyperfine fields in a quantitative manner and to estimate also the influence of next- 
nearest and further distant Cu atoms. The isotropic hyperfine constants at the oxygens 
were assumed to depend on the spin densities of the adjacent (NN), next-nearest-neighbors 
(NNN) and further distant Cu atoms according to 

D(0)= 7 E P(Cuj)+7' E P(Cuj)+7" E P( Cu i)- ( 9 ) 

jeNN jeNNN jgNNNN 
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For Cu, the isotropic hyperfine constants were decomposed into 

D(Cui) = ap(Cui) + /3 £ p(Cuj) +P' £ p(Cuj) + (3" £ p(Cuj). (10) 

jeNN jeNNN jeNNNN 

Figs. [8] and ^ already show that the main contributions are due to NN. This result is 
corroborated by least squares fits of all data to Eqs. (Q) and ([TO]) which yield 

a = -2.54 ±0.04 = 1.02 ± 0.02 | \ < 0.02 | (5" \ < 0.02 (11) 

and 

7 = 0.925 ± 0.006 \l'\< 0.006 | 7" | < 0.006 . (12) 

As expected, the contributions from Cu ions beyond NN are marginal. 

The analyses performed so far have demonstrated that all hyperfine coupling parameters 
are proportional to the sum of the atomic spin densities at the adjacent Cu ions. The 
essential point is that in this way also the values calculated for lower spin multiplicities 
could be accounted for. These contain complementary informations to the results of high- 
spin states which strongly improve the quality of the results as can be seen in Figs. H and 
|9]. Incorporating negative values of SjP(Cuj) into the evaluation renders a much better 
understanding of the various contributions. 

We are now in a position to extrapolate the data from the finite clusters to an extended 
system. Considering the three different values of p(Cui) in the largest cluster (see Table |V|), 
we expect p(Cu) in the homogeneous case to be 0.70. 

This leads to the following values: a = -1.778, b = 0.714, | V |< 0.014, c = 0.648, and 
\c'\< 0.004. 

It should be emphasized again, that the extracted values for the contact terms depend 
on the particular choice of the exchange- correlation functional. We have performed a few 
calculations with other functionals, too. They deliver values that differ by 15 % but they also 
lead to the conclusion that contributions from further distant neighbors are not important. 



B. Dipolar Hyperfine Interactions 

In contrast to the core-polarization, which is proportional to the spin density difference 
at a single point, the dipolar hyperfine coupling results from a spatial average of 1/r 3 with 
wave-functions. Its value, therefore, can be determined with much greater accuray than that 
of the contact interaction and does less crucially depend on the calculational method or the 
choice of the exchange-correlation functionals. It is given by 

4 P (R) = £ (V&tt I T ij (r - R) I V£(r)) - £ (t&(r) | T«(r - R) | ^,(r)) (13) 

m m ' 

where 



x 



3x;Xj — 5jjX 



2 



X 5 



(14) 
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In Table VII] the values for a dip (for Cu) and c dip (for 0) obtained with various clusters are 

given. It should be remarked that c dip for planar oxygen refers to the directions along the 
Cu-O-Cu bonds. The tensor coupling is only nearly axially symmetric. 

In analogy to the isotropic case, the dipolar contributions also turned out to depend 
mainly on the atomic spin densities of the NN copper ions besides the on-site term for Cu 
which is proportional to the on-site atomic spin density. A least squares fit of all data 
(including those with lower multiplicities) to the ansatz 

a di P ( Cu i) = "dipP(Cui) + /3 dip E P( Cu i) +/ 3 di P E p( Cu j) +#iip E p( Cu j) ( 15 ) 

jeNN jGNNN jGNNNN 

gives a dip = -5.206 ± 0.014, (3 dip = 0.101 ± 0.006, and | (3' dip |< 0.005. Thus there is also a 
small transferred dipolar hyperfine interaction for Cu. 
For the oxygens, the ansatz 



c. 



dip 



(O) = 7di P E P( Cu j) +7di P E P(Cuj)+7* P E P( Cu j)- ( 16 ) 



jeNN jGNNN jGNNNN 



leads to 7 dip = 0.527 ± 0.006, | 7 dip |< 0.006, and | 7 d ' ip |< 0.006. 

Performing the same extrapolation as was applied to the isotropic contributions, we 
finally get a dip = —3.644, b dip = 0.071, c^ = 0.369 (along the bond direction), c d f p = —0.177 
(perpendicular to the bond), and c d ; p = —0.192 (perpendicular to the plane). 

Represented in Table [IX] are the magnetic hyperfine couplings calculated for O(p). The 
spin-orbit contributions to the hyperfine fields are expected to be small in the case of O. 
We neglect them completely by assuming c" ot = c + c dip . These values mav be compared to 
those which have been extracted from various experiments in Refs.0 andS 



C. Origin of the Transferred Hyperfine Fields 

To investigate the mechanisms of spin transfer, we will analyze in detail the results from 
the two clusters containing one and two copper atoms and point out the differences. We 
commence by noting that the basic notions of spin-polarized density functional theoryl! rely 
on the concept of expressing all physical relevant quantities in terms of the spin densities 
alone. The "wave functions" are just auxiliary entities that are introduced to solve the 
Kohn-Sham equations^ and have no direct physical meaning. Nevertheless, to understand 
the origins of transferred hyperfine fields, it is convenient to assign the AOs which are built up 
from localized basis functions to atomic single-electron wave functions. Similarly, the MOs 
are interpreted as hybridizations between the individual AOs. Another remark concerns the 
representation of the AOs in terms of basis functions. We recall that the employed basis 
sets describe a "single-electron wave function" as a linear combination of several contracted 
radial Gauss functions. For Cu, there are three radial functions for each of the five 3d- 
orbitals and nine functions of s-like character. For O, we have four s-functions and three 
radial functions for each of the 2p-orbitals. 

Let us first consider the (CuO 6 ) 10 ~ ion embedded in the appropriate lattice. Overlap and 
covalent effects convey spin density from the Cu 2+ ion onto the ligand oxygen sites whose 
spin direction is parallel to that of the local Cu moment. The data in Table || show that in 
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the DF calculations the atomic spin density p(Cu) is reduced from 1 to 0.667 in favor of a 
spin density transfer to the four planar oxygens (4 x p(0(p)) = 0.328). The transfer from the 
Cu to the O(p) is much less pronounced at the HF level where p(Cu) = 0.902. This trend 
is also reflected in the dipolar hyperfme coupling (see Table |VII]| ): a^ p = 1.763 (2.210) 
and c£p = 0.388 (0.181) for DF (HF) with a total of ajjf p + 2 x c£ p = 2.539 (2.572). The 
transferred spin density is mainly on the O 2p CT orbital, but a small amount also goes to 
the O 2s which is thus expected to be polarized parallel to the spin density on the Cu. 
Indeed, we get D 2s (0) = 0.772, as can be seen in Table [X] where the core polarizations D(O) 
from DF calculations are listed for the individual orbitals. The 2s orbital in turn polarizes 
the Is but with opposite sign (D ls (0) = —.141). Since the MOs are linear combinations 
of AOs, the contact term may also have contributions that come from a product of an s- 
like AO centered at the nucleus under consideration with an AO centered on an adjacent 
atom. These contributions turn out to be small. They are listed separately as numbers in 
parentheses in Tables [X] and XI. On the Cu, the derealization of spin densities implies the 
following changes. The MO second highest in energy (at —.14 eV) is occupied with a spin- 
down electron. It mainly mixes the Cu 3d3 Z 2_ r 2 with the O(a) 2p z but has also some small 
contributions from the Cu 4s and the O(p) 2p a orbitals. The admixture of the 4s is larger 
than that of the following orbital (at —.58 eV) which has the same symmetry but is occupied 
with a spin-up electron. These two MOs thus cause the negative value D 4s (Cu) = —.650 (see 
Table [Xlf ). The 4s AO in turn polarizes the inner s orbitals. The alternating signs of the 4s, 
3s, and 2s contributions!! are a consequence of the Pauli exclusion principle. 

We next consider the system consisting of 2 Cu atoms with parallel moments and with 
1 bridging O (Ob r ) and 6 planar O (Oi) bonded to one of the Cu only. The DF calculations 
yield the atomic spin densities p(Cu) = 0.668, p(Ob r ) = 0.130, and p(Oi) = 0.089. Since the 
2p CT orbital can only share a maximal amount of hole contribution, the additional transfer 
of spin density onto the Ob r is less than twice the value in the CuOe system. The missing 
fraction is added to the six Oi. More positive spin density, however, can be put on the 
2s which leads to the total value D(Ob r ) = 1.275 (see Table [Xj). The dipolar contributions 
are a^f p = 1.746 and cgf p = 0.324. On the Cu atoms, the spin density distributions in the 
inner s shells cannot be changed significantly compared to the CuOg case since not much 
redistribution is possible in the inner core. Therefore the additional positive spin density 
resides mainly on the 4s atomic orbital. 

We thus get a coherent picture of the spin density distributions for a single Cu 2+ ion 
surrounded by four planar oxygens and of the additional spin transfers caused by a second 
copper ion. The interpretation in terms of localized AOs turned out to be helpful for the 
understanding of the spin densities which are the relevant quantities in density functional 
theory. The results of our quantum chemical calculations demonstrate that the spin densi- 
ties, when attributed to the Cu or to the O, are strongly connected. In the ground state, 
there is no reason for advocating a two-component model and we also see little chance that 
low-energy excited states could change this fact. 



D. Spin-Orbit Coupling 

For Cu, the spin-orbit coupling gives an appreciable contribution to the total hyperfine 
fields, especially in the perpendicular direction since accidentally a and a^ almost cancel 
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(a + a^ ip = 0.003). A calculation of the spin-orbit interaction cannot be carried out at the 
same level of quality as is possible for the other hyperfine interactions. Therefore we adopt 
the values for the isolated Cu 2+ ion from Ref.0 in the simplified form used in Ref.i: 

62 

k<r- 3 > (17) 



7 



and 



-yk<r~ 3 > (18) 



with a parameter value k = —0.044. With < r -3 >= 6.171, we get a| Q = 2.405 and a^, = 
0.427. It should be noted that the uncertainty in the spin-orbit coupling parameter and the 
excitation energies is comparatively large. Therefore, in Ref.S a range of values for k of ± 
20 % was considered. 



In Table [XI I| all our calculated values for the hyperfine coupling parameters for Cu are 
compiled. The spin-orbit contributions are put in parentheses to emphasize that they are 
estimated. 

No direct experimental information is available on these parameters. They can be deter- 
mined indirectly from a combination of anisotropic Knight shifts and relaxation rates and 
the nuclear resonance frequencies in antiferromagnetic compounds. The results of various 
analyses are shown in Table [XIl] . Our values are in reasonable agreement with the data. 



VI. SUMMARY AND CONCLUSIONS 

The electronic structure of La 2 Cu04 has been investigated by first-principles cluster 
calculations. These have the advantage that local properties can be studied in great detail 
but usually suffer from a somewhat uncontrolled embedding in the periodic lattice and 
surface problems. By using a sequence of clusters containing up to nine copper atoms, we 
have demonstrated which quantities can be evaluated with clusters of modest size. 

As concerns the theoretical approximations to treat the many-electron problem in Cu02 
planes, we have shown that distinct differences between the Hartree-Fock and density func- 
tional methods exist. From our results we conclude that the former approach, which neglects 
correlations entirely, is not adequate for cuprate superconducting materials. 

We have evaluated the electric-field gradients at the Cu and the O sites for a variety of 
cluster sizes and spin multiplicities. For both planar and apical oxygens a satisfactory agree- 
ment with the experimentally determined quadrupole frequencies was obtained. For copper, 
a comparison between theoretical and experimental values is hampered by the uncertainty 
with which the Cu nuclear quadrupole moments are known. 

With respect to magnetic hyperfine properties one has to distinguish between the core- 
polarization and the dipolar contributions. The latter can be determined theoretically with 
good reliability. The former, however, is a quantity whose value is given by subtle cancel- 
lations of contributions from various s-like atomic orbitals. Our calculations demonstrated 
that besides the negative on-site core-polarization for Cu, a sizeable positive contribution 
transferred from neighboring Cu + ions exists as has been suggested by Mila and Rice!. 
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We have investigated this transfer of spin density in detail and found that there is no 
appreciable contribution from copper ions other than the nearest neighbors. This is the main 
result of the present work since it has far reaching consequences. In particular, it questions 
the reconciliation between NMR and neutron scattering experiments and points out that the 
disagreements are now as before. By introducing a transferred hyperfine coupling between 
next-nearest neighbor Cu spins and 17 nuclei spin, Zha et alJliJ were able to conciliate 
many different experiments. We have shown that there is no microscopic justification for 
the presence of these additional hyperfine couplings. The discrepancies among the data 
therefore need another explanation. 

On the other hand, the results of the analysis of the origins of the calculated hyperfine 
fields carried out in Sec. |VQ shows that the spin densities on the coppers and oxygens are 
very tightly connected. This gives strong support for a one-component model of the spin 
fluid, at least as concerns ground-state properties. 

We have reached these conclusions by carrying out extended ab-initio cluster calculations. 
One has, of course, to take into consideration that this approach to describe the electronic 
structure emphasizes local features and is surely far more appropriate for insulating materials 
than for metals. It would be a big surprise, however, when a tight-binding picture completely 
failed to reproduce qualitative aspects of a system with itinerant charge carriers. This would 
shift the problems to a quite different field. 

The calculated values for the magnetic hyperfine couplings are in agreement with those 
extracted from various NMR and NQR experiments if we adopt the conventional estimates 
for the spin-orbit interaction for the Cu 2+ ion. In this respect, improved quantum chemical 
calculations of the spin-orbit coupling would be desirable. Although our results give a small 
value for the hyperfine field that contributes to the spin part of the Knight shift in c-direction 
(a" + 4b + 4b"), we think that the vanishing of Kjj below T c in all high-temperature super- 
conductors and with all doping concentrations is still another open problem that deserves 
further experimental and theoretical consideration. 
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FIGURES 



FIG. 1. The CuOe/Cu4Laio cluster. 

FIG. 2. Energies of the highest occupied MOs in the CuOe/C^Laio cluster with contributions 
from the individual AOs: (a) HF, (b) DF. The length of the bar is proportional to the square of 
the expansion coefficient of the MO into the corresponding AOs. Spin-up (spin-down) orbitals are 
denoted by solid (dashed) bars. 

FIG. 3. Total and partial density of states of the highest "band" in the CuOe/Cu4Laio cluster 
as obtained with the DF method. 

FIG. 4. The CugC^/Cu^Laso cluster. 

FIG. 5. Highest occupied molecular orbital for the Cug042/Cui2La5o cluster. 

FIG. 6. Difference between spin-up and -down densities at Cu vs. number of nearest neighbor 
Cu atoms, N, for different sites in various clusters with n Cu atoms. Circles: n = 1, diamonds: 
n = 3, triangle down: n = 4, triangle up: n = 5, squares: n = 9. DF (HF) values correspond to 
the upper (lower) data set. 

FIG. 7. Difference in spin densities along the O-Cu-O-Cu-O-Cu-0 bonds as calculated with 
DF/GGA for the cluster Cus026/Cu8La34. From top to bottom: multiplicity m = 6, m = 4, and 
m = 2 horizontally (vertically) with spin density signs H h ( h). 

FIG. 8. Difference between spin-up and -down densities at O(p) vs. sum of atomic spin 
densities at neigboring Cu atoms for different sites in various clusters: Cug (diamonds), Cus, m=6 
and m=4 (circles), CU4 (squares), and Cus with m=2 (open circles). 

FIG. 9. D(Cui)/p(Cui) values plotted against the sum of atomic spin densities at the NN Cu 
atoms divided by p(Cui): Cug (diamonds), Cus, m=6 and m=4 (circles), and Cus with m=2 (open 
circles) . 
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TABLES 



TABLE I. Clusters, number of atoms N, atoms with pseudopotentials PP, total number of 
electrons E, number of basis functions B, and number of primitive Gaussian functions P. 



cluster N PP E B P 

Cu0 6 7 Cu 4 Lai 87 117 222 

Cu 2 0n 13 Cu 6 Lai2 164 221 418 

Cu 4 O 20 24 Cu 8 La 26 308 334 772 

Cu 5 26 31 Cu 8 La 34 395 533 1006 

CugO-42 51 Cui 2 La 50 663 621 1116 



TABLE II. EFG component V Z z, core polarization a cp , dipolar contribution to the hyperfine 
tensor a|j , and atomic spin density p at the Cu, as obtained with HF and DF/GGA calculations 
for the cluster CuOe/Ch^Laio- 



method 


basis 


v zg 






P 


HF a 




2.232 


-0.503 


-4.421 


.88 


HF 6 




1.840 








HF 


6-311G 


1.936 


-3.519 


-4.421 


.902 


DF/BLYP 


6-311G 


1.419 


-1.784 


-3.526 


.667 


DF/BLYP 


6-311Gpd 


1.360 


-1.784 


-3.502 


.661 


DF/SVWN 


6-311G 


1.364 


-1.625 


-3.425 


.642 


DF/XALYP 


6-311G 


1.367 


-1.801 


-3.409 


.635 


DF/BPW91 


6-311G 


1.415 


-1.943 


-3.517 


.669 



a Ref . [2g 
b Ref. 
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TABLE III. Contributions to the EFG component V zz , calculated with DF/GGA for the 
cluster CuOg/Cu4Laio, from the nuclei within the cluster, from the point charges around the 
cluster, and from the individual shells. The "remainder" lists just the (small) rest of all shells not 
explicitly given in the table. 



Nuclei 






0.375 


Point charges 






0.012 




spin up 


spin down 


sum 


Vxi Vy 


-1.183 


-1.210 


-2.393 


Vz 


0.384 


0.386 


0.770 


d-x 2 — y 2 i d X y 


-9.443 


-5.720 


-15.163 


d 3z 2_ r 2 


4.446 


4.381 


8.827 


d-xzi dyz 


4.522 


4.481 


9.003 


Remainder 


0.027 


-0.039 


-0.012 


Total 






1.419 



TABLE IV. EFG component \ zz and asymmetry parameter r/ for the Cu site in the La2Cu04 
system calculated for different cluster sizes and various spin multiplicities m. The basis set used 
for all calculations was 6-311G. For comparison, the theoretical results of Sulaiman et al.@ and 
Martinll, as well as the measured quadrupole frequency^, have been transformed into atomic 

units. 

HF LDA/GGA 



cluster 


m 


v 2 , 


V 


Vzz 


V 


CuOe/Cu4Laio 


2 


1.936 





1.419 





Cu 2 0ii/Cu 6 Lai2 


3 


1.904 


0.04 


1.360 


0.14 


Cu 2 0ii/Cu 6 Lai2 


1 


1.884 


0.04 


1.181 


0.14 


Cu 4 02o/Cu 8 La 2 6 


5 


1.925 


0.00 


1.354 


0.01 


Cu 5 026/Cu 8 La34 


6 


1.963 





1.128 





Cu 5 026/Cu 8 La 3 4 


4 






1.071 





Cu 5 026/Cu 8 La34 


2 






1.097 


0.28 


Cu 9 O42/Cui2La40 


10 


1.975 





1.264 





CuOg 




2.234 


0.02 






CuOe/CuiLaJo 




1.840 








Cu 2 0ii/Cu 6 La 6 12 




1.816 








experiment 




1.560 ( 63 Q 


= -0.18 b) 


1.331 


( 63 Q = -0.211 b) 



a Ref . |3 
b Ref. £2 
c Ref. Wi 
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TABLE V. EFG components 


tor the ptanar oxy 


gen 


0(p). 








HF 






LDA/GGA 




cluster 


V 


^yy 


v zz 






Vzz 


Cu 2 0n/Cu 6 Lai2 


-.416 


.335 


.081 


.858 


.518 


.340 


Cu 4 02o/Cu 8 La 2 6 


-.403 


.363 


.040 


.839 


.536 


.303 


Cu 5 026/Cu 8 La 3 4 


-.411 


.326 


.049 


.847 


.545 


.302 


Cu 9 042/Cui2La 5 o 


-.499 


.353 


.146 


.881 


.554 


.327 


experiment" 


-.75 


.51 


.24 


-.75 


.51 


.24 



Ref.0 



TABLE VI. Atomic spin densities and D-values at different planar sites in the cluster 
Cug042/Cui2La5o as obtained with the DF method. N denotes the number of nearest neighbor 

Cu ions. 

position N p D 

Cu(0/0) 4 .6935 0.9215 

Cu(2/0) 3 .6870 0.3209 

Cu(2/2) 2 .6766 -0.3636 

O(l/0) 2 .1427 1.2550 

0(2/1) 2 .1390 1.2214 

O(3/0) 1 .0896 0.6501 

0(3/2) 1 .0815 0.6384 



TABLE VII. Atomic spin densities and D-values at different planar sites in the cluster 
Cu5026/CugLa34 as obtained for spin multiplicities m = 6,4, and 2 with the DF method. 



m=6 m=4 m=2 



position 


P 


D 


P 


D 


P 


D 


Cu(0/0) 


.6745 


1.1963 


-.2636 


3.2337 


-.3036 


1.5959 


Cu(2/0) 


.6743 


-0.9500 


.5738 


-1.7329 


.5412 


-1.6747 


O(l/0) 


.1356 


1.2034 


.0260 


0.3062 


.0199 


0.2086 


O(3/0) 


.0818 


0.6258 


.0726 


0.5500 


.0677 


0.5228 


0(2/1) 


.0872 


0.6417 


.0705 


0.5603 


.0660 


0.4520 


Cu(0/2) 










-.6364 


1.3237 


Cu(0/-2) 










.5310 


-1.6093 


O(0/l) 










-.0893 


-0.8956 


O(0/3) 










-.0778 


-0.6023 


O(0/-l) 










.0143 


0.3083 


O(0/-3) 










.0697 


0.5043 


0(1/2) 










-.0773 


-0.5994 


0(2/-l) 










-.0651 


0.5345 


0(l/-2) 










.0657 


0.5265 
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TABLE VIII. Values for the dipolar hyperfine couplings for Cu and 0(p) for different clusters. 



cluster a^ p (Cu) c d i P ( Q ) 

HF DF HF DF 

CuO 6 /Cu 4 Lai -4.421 -3.526 .181 .388 

Cu 2 0n/Cu 6 Lai2 -4.430 -3.493 .322 .648 

Cu 4 2 o/Cu 8 La26 -4.452 -3.479 .322 .656 

Cu 5 026/Cu 8 La 3 4 -4.486 -3.249 .324 .668 

Cu 9 042/Cui2La 4 o -4.491 -3.370 .324 .684 



TABLE IX. Compilation of results for the hyperfine interactions at O(p) and comparison with 
values extracted from experiments. 





c 


P xx 
c dip 




„zz 


„xx 

Hot 


r yy 
Hot 


„zz 

Hot 


present 


0.648 


0.369 


-0.177 


-0.192 


1.017 


0.471 


0.456 


La 2 Cu0 4 a 










0.748 


0.480 


0.527 


La 2 Cu0 4 b 










0.863 


0.511 


0.615 



a Ref. y 
b Ref. Wl 



TABLE X. Contributions to the core-polarization at an oxygen, D(O), with one (N=l) adjacent 
Cu ion, between two (N=2) adjacent Cu ions, and difference. The numbers in parentheses are 



contributions from AO 


s centered 


at 


neighboring 


nuclei. 








atomic orbital 


N = 


1 




N = 2 




Difference 




Is 


-0.141 




(.000) 


-0.229 


(.001) 


-0.088 ( 


.001) 


2s 


0.772 




(.012) 


1.503 


(.052) 


0.731 ( 


.041) 



Total 0.632 (.012) 1.274 (.053) 0.643 (.042) 



TABLE XL Contributions to the core-polarization at a copper, D(Cu), with no adjacent Cu 
ions (N=0), with one NN Cu ion, and their difference. The numbers in parentheses are contributions 
from AOs centered at neighboring nuclei. 



atomic orbital 


N = 





N = 1 




Difference 




Is 


-0.060 


(.000) 


-0.059 


(.000) 


0.001 


^.000) 


2s 


-3.733 


(.001) 


-3.687 


(.002) 


0.046 


(.001) 


3s 


2.659 


(.009) 


2.618 


(.009) 


-0.041 


^.000) 


4s 


-0.650 


(.001) 


-0.029 ( 


-.003) 


0.621 (- 


-.004) 


Total 


-1.784 


(.011) 


-1.157 


(.008) 


0.627 (- 


-.003) 
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TABLE XII. Compilation of results for the hyperfine interactions at Cu and comparison with 
values extracted from experiments. 



b 



a dip ^dip ^f°_ 



a 11 

"tot 


a tot 


3.017 


0.471 


2.955 


0.288 


1.550 


0.447 


3.636 


-0.075 


3.115 


0.400 


3.010 


-0.179 


2.748 


0.495 


2.524 


0.591 



present 
La 2 Cu0 4 a 
La 2 Cu0 4 b 
YBa 2 Cu 3 7 c 
YBa 2 Cu 3 7 d 
YBa 2 Cu 3 7 e 
YBa 2 Cu 3 7 a 
YBa 2 Cu 3 7 b 

a Ref. [11 
b Ref. |6 
c Ref. | 
d Ref. | 
e Ref. | 



-1.778 



-2.050 



0.714 
0.577 
0.607 
0.744 
0.585 
0.712 
0.687 
0.623 



-3.644 0.071 (2.405) (0.427) 



-3.672 



2.086 0.289 
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9 • 9 



E-E F 
[eV] 


Cu 


0(p) 


0(a) 


dzz d X z d vz d xx d xv 


Px Pv Pz 
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Cu 


0(p) 


0(a) 


[eV] 
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Px Pv Pz 


Px Pv Pz 


2 - 
1 - 
- 
-1 - 
-2 - 














-3 - 
-4 - 




m m 




-5 - 
-6 - 
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